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Abstract 

We consider a class of nonlinear Schrodinger / Gross-Pitaveskii (NLS-GP) equa- 
tions, i.e. NLS with a linear potential. We obtain conditions for a symmetry breaking 
bifurcation in a symmetric family of states as M, the squared norm (particle number, 
optical power), is increased. In the special case where the linear potential is a double- 
well with well separation L, we estimate J\fcr{L), the symmetry breaking threshold. 
Along the "lowest energy" symmetric branch, there is an exchange of stability from 
the symmetric to asymmetric branch as M is increased beyond Mcr- 
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1 Introduction 



Symmetry breaking is a ubiquitous and important phenomenon which arises in a wide range 
of physical systems. In this paper, we consider a class PDEs, which are invariant under a 
symmetry group. For sufficiently small values of a parameter, A/", the preferred (dynami- 
cally stable) stationary (bound) state of the system is invariant under this symmetry group. 
However, above a critical parameter, Mcri although the group-invariant state persists, the 
preferred state of the system is a state which (i) exists only for M > Mct and (ii) is no longer 
invariant. That is, symmetry is broken and there is an exchange of stability. 

Physical examples of symmetry breaking include liquid crystals [23], quantum dots |27j . 
semiconductor lasers [S] and pattern dynamics [23]. This article focuses on spontaneous sym- 
metry breaking, as a phenomenon in nonlinear optics [SlilSlill], as well as in the macroscopic 
quantum setting of Bose-Einstein condensation (BEC) [T] . Here, the governing equations are 
partial differential equations (PDEs) of nonlinear Schrodinger / Gross-Pitaevskii type (NLS- 
GP). Symmetry breaking has been observed experimentally in optics for two-component 
spatial optical vector solitons (i.e., for self-guided laser beams in Kerr media and focusing 
cubic nonlinearities) in [3], as well as for the electric field distribution between two- wells 
of a photorefractive crystal in [16] (and between three such wells in [H]). In BECs, an 
effective double well formed by a combined (parabolic) magnetic trapping and a (periodic) 
optical trapping of the atoms may have similar effects [1] , and lead to "macroscopic quantum 
self-trapping" . 

Symmetry breaking in ground states of the three-dimensional NLS-GP equation, with 
an attractive nonlinearity of Hartree-type and a symmetric double well linear potential, was 
considered in Aschbacher et. al. [2]; see also Remark 12. 1[ Ground states are positive and 
symmetric nonlinear bound states, arising as minimizers of, 7Y, the NLS-GP Hamiltonian 
energy subject to fixed. A/", the squared norm. For the class of equations considered in [2], 
ground states exist for any A/" > 0. It is proved that for sufficiently large A/", any ground state 
is concentrated in only one of the wells, i.e. symmetry is broken. The analysis in [2] is an 
asymptotic study for large A/", showing that if M is sufficiently large, then it is energetically 
preferable for the ground state to localize in a single well. In contrast, at small norm 
the ground state is bi-modal, having the symmetries of the linear Schrodinger operator 
with symmetric double-well potential. For macroscopic quantum systems, the squared LP' 
norm, denoted by A/", is the particle number, while in optics it is the optical power. An 
attractive nonlinearity corresponds to the case of negative scattering length in BEC and 
positive attractive Kerr nonlinearity in optics. 

An alternative approach to symmetry breaking in NLS-GP is via bifurcation theory. It 
follows from [211 EO] that a family of "nonlinear ground states" bifurcates from the zero 
solution {M = 0) at the ground state energy of the Schrodinger operator with a linear 
double well potential. This nonlinear ground state branch consists of states having the 
same bi-modal symmetry of the linear ground state. In this article we prove, under suitable 
conditions, that there is a secondary bifurcation to an asymmetric state at critical A/" = 
A/'cr > 0. Moreover, we show that there is a transfer or exchange of stability which takes 
place at A/'cr-; for A/" < N'cr the symmetric state is stable, while for A/" > A/'cr the asymmetric 
state is stable. Since our method is based on local bifurcation analysis we do not require 
that the states we consider satisfy a minimization principle, as in [2] . Thus, quite generally. 



2 



symmetry-breaking occurs as a consequence of the (finite dimensional) normal form, arising 
in systems with certain symmetry properties. Although we can treat a large class of problems 
for which there is no minimization principle, our analysis, at present, is restricted to small 
norm. As we shall see, this can be ensured, for example, by taking the distance between 
wells in the double-well, to be sufficiently large. 

In [To] the precise transition point to symmetry breaking, A/'cr, of the ground state and 
the transfer of its stabihty to an asymmetric ground state was considered (by geometric 
dynamical systems methods) in the exactly solvable NLS-GP, with a double well potential 
consisting of two Dirac delta functions, separated by a distance L. Additionally, the behavior 
of the function Afcr{L), was considered. Another solvable model was examined by numerical 
means in [18]. A study of dynamics for nonlinear double wells appeared in |22j . 

We study Afcr{L), in general. Afcr{L), the value at which symmetry breaking occurs, is 
closely related to the spectral properties of the linearization of NLS-GP about the symmetric 
branch. Indeed, so long as the linearization of NLS-GP at the symmetric state has no non- 
symmetric null space, the symmetric state is locally unique, by the implicit function theorem 
[T9] . The mechanism for symmetry breaking is the first appearance of an anti-symmetric 
element in the null space of the linearization for some A/" = A^r- This is demonstrated for a 
finite dimensional Galerkin approximation of NLS-GP in [THIIIB]- The present work extends 
and generalizes this analysis to the full infinite dimensional problem using the Lyapunov- 
Schmidt method [19]. Control of the corrections to the finite-dimensional approximation 
requires small norm of the states considered. Since, as anticipated by the Galerkin approx- 
imation, N'cr is proportional to the distance between the lowest eigenvalues of the double 
well, which is exponentially small in L, our results apply to double wells with separation L, 
hold for L sufficiently large. 

The article is organized as follows. In section [2] we introduce the NLS-GP model and 
give a technical formulation of the bifurcation problem. In section [3] we study a finite 
dimensional truncation of the bifurcation problem, identifying a relevant bifurcation point. 
In section HI we prove the persistence of this symmetry breaking bifurcation in the full 
NLS-GP problem, for A/" > N'cr- Moreover, we show that the lowest energy symmetric 
state becomes dynamically unstable at A^r and the bifurcating asymmetric state is the 
dynamically stable ground state for A/" > J^cr- Figure [1] shows a typical bifurcation diagram 
demonstrating symmetry breaking for the NLS-GP system with a double well potential. At 
the bifurcation point N'cr (marked by a circle in the figure), the symmetric ground state 
becomes unstable and a stable asymmetric state emanates from it. 

The main results are stated in Theorem 14.11 Corollary 14.11 and Theorem 15.11 In partic- 
ular, we obtain an asymptotic formula for the critical particle number (optical power) for 
symmetry breaking in NLS-GP, 



Here, (f2o,^/'o) a-^d (Qijtpi) are eigenvalue - eigenfunction pairs of the linear Schrodinger 
operator H = —A + V, where Qq and Qi are separated from other spectrum, and S is 
a positive constant, given by (14.11) . depending on ipQ and ipi. The most important case is 
where Qq < Qi are the lowest two energies (linear ground and first excited states). For 
double wells with separation L, we have A/'cr = NcriL), depending on the eigenvalue spacing 




(1.1) 
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Figure 1: (Color Online) Bifurcation diagram for NLS-GP with double well potential (16.11) 
with parameters s = 1, L = Q and cubic nonlinearity. The first bifurcation is from the 
the zero state at the ground state energy of the double well. Secondary bifurcation to an 
asymmetric state at A/" = A/'cr is marked by a (red) circle. For A/" < N'cr the symmetric state 
(thick (blue) solid line) is nonlinearly dynamically stable. For A/" > N'cr the symmetric state 
is unstable (thick (blue) dashed line). The stable asymmetric state, appearing for M > Mcr, 
is marked by a thin (red) solid line. The (unstable) antisymmetric state is marked by a thin 
(green) dashed line. 

^q{L) — r2i(L), which is exponentially small if L is large and S is of order one. Thus, for 
large L, the bifurcation occurs at small norm. This is the weakly nonlinear regime in 
which the corrections to the finite dimensional model can be controlled perturbatively. A 
local bifurcation diagram of this type will occur for any simple even-odd symmetric pair of 
simple eigenvalues of H in the weakly nonlinear regime, so long as the eigen-frequencies are 
separated from the rest of the spectrum of if; see Proposition 14. II and the Gap Condition 
(14. 7p . Therefore, a similar phenomenon occurs for higher order, nearly degenerate pairs of 
eigen-states of the double wells, arising from isolated single wells with multiple eigenstates. 
Section [6] contains numerical results validating our theoretical analysis. 
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DMS-0405921 and DMS-060372. PGK was supported, in part, by DMS-0204585, NSF- 
CAREER and DMS-0505663, and acknowledges valuable discussions with T. Kapitula and 
Z. Chen. MIW was supported, in part, by DMS-0412305 and DMS-0530853. Part of this 
research was done while Eli Shlizerman was a visiting graduate student in the Department 
of Applied Physics and Applied Mathematics at Columbia University. 
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2 Technical formulation 



Consider the initial-value-problem for the time-dependent nonlinear Schrodinger / Gross- 
Pitaevskii equation (NLS-GP) 

idt'ip = Hip + g{x)K[ipip ] ^/^(x, 0) given (2.1) 
H = -A + V{x). (2.2) 

We assume: 

(HI) The initial value problem for NLS-GP is well-posed in the space C°([0, oo); if^(lR")). 

(H2) The potential, V{x) is assumed to be real-valued , smooth and rapidly decaying as 
\x\ — >■ oo. The basic example of V{x), we have in mind is a double- well potential, consisting 
of two identical potential wells, separated by a distance L. Thus, we also assume symmetry 
with respect to the hyperplane, which without loss of generality can be taken to be {xi = 0}: 

V{Xi,X2,...,Xn) = V{-Xi,X2,...,Xn). (2.3) 

We assume the nonlinear term, K[iljip], to be attractive, cubic ( local or nonlocal), and 
symmetric in one variable. Specifically, we assume the following 
(H3) Hypotheses on the nonlinear term: 

(a) g{xi, X2,..., Xn) = g{-xi, X2, • • • , Xn) (symmetry) 

(b) g{x) < (attractive / focusing) 

(c) K[h] = J K{x - y)h{y)dy, K{xi, X2, . . . , x^) = K{-xi, X2, . . . , x^), K > 0. 

(d) Consider the map N : H"^ x H"^ x H"^ ^ L"^ defined by 

iV(01,02,03)=^7/^[0102]03. (2.4) 

We also write N{u) = N{u,u,u) and note that duN{u) = N{-,u,u) + N{u,-,u) + 
N{u,u, ■). We assume there exists a constant k > such that 

||iV(0O,01,02)||L2 < k \\Mh4<P2\\h4<P3\\h^- (2.5) 

Several illustrative and important examples are now given: 

Example 1: Gross-Pitaevskii equation for BECs with negative scattering 

length g{x) = —1, K{x) = 6{x) 

Example 2: Nonlinear Schrodinger equation for optical media with a nonlocal kernel 
g{x) = ±1, K{x) = Aexp(— x^/o"^) [I7j (see also |1] for similar considerations in BECs). 
Example 3: Photorefractive nonlinearities The approach of the current paper can be 
adapted to the setting of photorefractive crystals with saturable nonlinearities and appro- 
priate optically induced potentials [5]. The relevant symmetry breaking phenomenology is 
experimentally observable, as shown in [16]. 
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Nonlinear bound states: Nonlinear bound states are solutions of NLS-GP of the form 

= e-^™*n(x), (2.6) 

where E F^(1R") solves 

H^n + 9{x) K[\^n\''] - fi^n = 0, uE H' (2.7) 

If the potential V{x) is such that the operator H = —A + V{x) has a discrete eigenvalue, 
-E*, and correspsonding eigenstate ijj^, then for energies near E near E^, and one expects 
small amplitude nonlinear bound states, which are to leading order small multiples of 
This is the standard setting of bifurcation from a simple eigenvalue [19] , which follows from 
the implicit function theorem. 

Theorem 2.1 fWj l21 f Let (\1/, E) = (■?/'*, E^,) be a simple eigenpair, of the eigenvalue problem 
H'^ = i.e. dim{p : {H — E^)p = 0} = 1. Then, there exists a unique smooth curve of 
nontrivial solutions a ^ ( a),Q{a) ), defined in a neighborhood ofa = 0, such that 

= a{^o + Oi\a\^)), = fio + Oi\a\^), a ^ 0. (2.8) 

Remark 2.1 For a large class of problems, a nonlinear ground state can be characterized 
variationally as a constrained minimum of the NLS / GP energy subject to fixed squared 
norm. Define the NLS-GP Hamiltonian energy functional 

UNLs^GPm = J iv$r + ^i$r dy + I j MJ^m'] dy (2.9) 

and the particle number (optical power) 

A/'[$] = j |$|2 dy, (2.10) 

where 

ICm^] = j K{x-y)\^{x)\^\^{y)\^ dy. (2.11) 
In particular, the following can be proved: 

Theorem 2.2 Let I\ = mij^[f]=x '^[/]- // ~oo < Ix < 0, then the minimum is attained 
at a positive solution of 1^2. 7| j. Here, Vt = f2(A) is a Lagrange multiplier for the constrained 
variational problem. 

In m the nonlinear Hartree equation is studied; K,[h] = \y\~^ -kh, g = —1. It is proved that if 
V{x) is a double-well potential, then for A sufficiently large, the minimizer does not have the 
same symmetry as the linear ground state. By uniqueness, ensured by the implicit function 
theorem, for small M , the minimizer has the same symmetry as that as the linear ground 
state and has the expansion ^2. 9) : see and section\^ 

We make the following 
Spectral assumptions on H 
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(H4) H has a pair of simple eigenvalues VLq and VLi. ipQ and ipi, the corresponding (real- 
valued) eigenf unctions are, respectively, even and odd in Xi. 



Example 2.1 The basic example: Double well potentials 

A class of examples of great interest is that of double well potentials. The simplest example, 
in one space dimension, is obtained as follows; see section\^for the multidimensional case. 
Start with a single potential well (rapidly decaying as \x\ ^ oo), Vo{x), having exactly one 
eigenvalue, u, Hqi/j^ = (—A + Vo{x))'ip^ = ivip^; see Figure ^^). Center this well at x = —L 
and place an identical well, centered at x = L. Denote by Vl{x) the resulting double-well 
potential and denote the Schrddinger operator: 



Hl = -A + Vl{x) 



(2.12) 



There exists L > Lq, such that for L > Lq, Hi has a pair of eigenvalues, VLq = Qo{L) and 
Qi = Qi{L), Qq < Qi, and corresponding eigenfunctions ipo and ipi; see Figure ^B>)- ipo is 
symmetric with respect to x = and ipi is antisymmetric with respect to x = 0. Moreover, 
for L sufficiently large, \flo — Qi\ = C(e~'*^), n > 0; see J^; see also section\^ 



(a) 



(b) 





-0.2 



-0.1 



0.1 



Figure 2: This figure demonstrates a single and a double well potential and the spectrum of 
H and Hl respectively. Panel (a) shows a single well potential and under it the spectrum of 
H, with an eigenvalue marked by a (red) mark 'o' at uj and continuous spectrum marked by 
a (black) line for energies u; > 0. Panel (b) shows the double well centered at ±L and the 
spectrum of underneath. The eigenvalues VLq and Vti are each marked by a (blue) mark 
'*' and a (green) mark 'x' respectively on either side of the location uj - (red) mark 'o'. The 
continuous spectrum is marked by a (black) line for energies > 0. 
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The construction can be generalized. If —A + Vq{x) has m bound states, then forming a 
double well Vl, with L sufficiently large, = —A + Vl will have m— pairs of eigenvalues: 
(f22j; ^2j+i)i J = 0,...,m — 1, eigenfunctions ip2j (symmetric) and ip2j+i anti- symmetric. 



By Theorem I2.H for small Af, there exists a unique non-trivial nonlinear bound state, 
bifurcating from the zero solution at the ground state energy, Qq, of H- By uniqueness, 
ensured by the implicit function theorem, these small amplitude nonlinear bound states 
have the same symmetries as the double well; they are bi-modal. We also know from [2] that 
for sufficiently large Af the ground state has broken symmetry. We now seek to elucidate 
the transition from the regime of Af small to A/" large. 

We work in the general setting of hypotheses (H1)-(H4). Define spectral projections 
onto the bound and continuous spectral parts of H: 

Po=(^o,-)^o, Pi = (V'i,-)V'i, P = I - Po - Pi (2.13) 

Here, 

{f,9) = I fgdx. (2.14) 



We decompose the solutions of Eq. (12. 7p according to 

= coipo + CiV^i + r], r] = Pr]. (2.15) 

We next substitute the expression (I2.15P into equation (12. 7p and then act with projections 
Pq, Pi and P to the resulting equation. Using the symmetry and anti-symmetry properties 
of the eigenstates, we obtain three equations which are equivalent to the PDE (12.71) : 

[Qq -Q)co + aoooolcopco + (aono + aooii) |cipCo + aooiiCiCo + {ipog,7l{co, Ci,r])) = 

(2.16) 

{Qi -Q)ci + aiiiilcipci + (aioio + aiooi) |copCi + aiowclci + {ipig,Tl{cQ,ci,ri)) = 

(2.17) 

(H-n) f] = -P g[F{-,co,ci) + 7^(co,Cl,r/) ] (2.18) 

F(-, Co, Ci) is independent of r] and 7^(co, Ci, rj) contains linear, quadratic and cubic terms in 
r]. The coefficients akimn are defined by: 

{^k,gK[^i^J^^n) (2.19) 

We shall study the character of the set of solutions of the system (12.161) . (I2.17p . (12.181) 
restricted to the level set 

\^u\^dx = M ^ |coP + |cip + / \r]\'^dx = M (2.20) 



as Af varies. 

Let VLq and VLi denote the two lowest eigenvalues of Hi. We prove (Theorem 14. 11 Corollary 
liTl Theorem EJ]): 
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• There exist two solution branches, parametrized by A/", which bifurcate from the zero 
solution at the eigenvalues, Qq and Qi. 

• Along the branch, {Q, emanating from the solution {Q = Qq, \1/ = 0) , there is a 
symmetry breaking bifurcation at A/" = J^crit > 0. In particular, let Ucrit denote the 
solution of (12. 7p corresponding to the value A/" = Merit- Then, in a neighborhood Ucrit, 
for Af < Merit there is only one solution of (12.71) . the symmetric ground state, while for 
M > Merit there are two solutions one symmetric and a second asymmetric. 

• Exchange of stability at the bifurcation point: For M < Merit the symmetric state 
is dynamically stable, while for M > Merit the asymmetric state is stable and the 
symmetric state is exponentially unstable. 



3 Bifurcations in a finite dimensional approximation 

It is illustrative to consider the finite dimensional approximation to the system fl2.16|2.17|2.18p . 
obtained by neglecting the continuous spectral part, Pu. Let's first set ?7 = 0, and therefore 
TZ{co, Ci, 0) = 0. Under this assumption of no coupling to the continuous spectral part of H, 
we obtain the finite dimensional system: 

(1^0 — f2) Co + aoooo|copco + (aoiio + Oooii) |cipCo + aooiic^co = (3-1) 
(fii - fi) ci + aiiii|cipci + (aioio + aiooi) |copci + aioioCoCi = (3.2) 
|coP + |cip = M (3.3) 

Our strategy is to first analyze the bifurcation problem for this approximate finite-dimensional 
system of algebraic equations. We then treat the corrections, coming from coupling to the 
continuous spectral part of H, rj, perturbatively. 

For simplicity we take Cj real: Cj = pj G IR; see section HI Then, 

po [r^o-^ + aooooPo + (aoiio + 2aooii)Pi ] = (3.4) 
pi + aiiiiPi + (fliooi + 2aioio) Po ] = (3.5) 

pI + pI - M = 0. (3.6) 



Introduce the notation 
Then, 



Vo = pI V^ = pI (3.7) 



J^o{Vo,Vi,n;M) = Vo [no-n + aoooo^^o + (aoiio + 2aooii)Pi ] = 

Ti{Vo,Vin-M) = Vi [^i-n + aiiiiPi + (aiooi + 2aioio) T'o ] = 
TM{Vo,Vr,n-M) =Vo + Vi - M = 0. {3.i 

Solutions of the approximate system 
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(1) Q^^^Af) = ( ) = ( AT, 0,^0 + aooooAT ) - approximate nonlinear 
ground state branch 

(2) Q(i)(7V) = ( ) = ( 0,7Vr,fii + aiiiiAr ) - approximate nonlinear 
excited state branch 

Thus we have a system of equations J-'{Q,Af) = 0, where JF : (1R^„ x 1R_|_ x IR) x 1R+ 
lR^xlR+, mapping (Q,7V) ^ J^(Q, AT) smoothly. We have that J^(Q(j)(7V),7V) = 0, j = 0, 1 
for all A/" > 0. A bifurcation (onset of multiple solutions) can occur only at a value of AC 
for which the Jacobian (iQjF(Q(-')(A/'*); AC) is singular. The point (Q(-'^(A/'*); AC) is called a 
bifurcation point. In a neighborhood of a bifurcation point there is a multiplicity of solutions 
(non- uniqueness) for a given A/". The detailed character of the bifurcation is suggested by 
the nature of the null space of (igjF(Q(-^)(AC); AC)- 

We next compute dQj-'{Q^^\Af);JV) along the different branches in order to see whether 
and where there are bifurcations. 



The Jacobian is given by 

r^o — + 2aoooo^o + (^0110 + aooii)'Pi (aoiio + 2aooii)^o —Po\ 

(aiooi + 2aioio)'Pi Qi - Q + 2aiiiiPi + (oono + 2aioio)'Po -'Pi (3.9) 

1 1 / 

A candidate value of A/" for which there is a bifurcation point along the "ground state 
branch" is one for which 

det(de^(Q(°)(Ar);Ar)) =0 ^ AT = A/:(o) ^ ^^1^%^ (3.10) 

c^oooo ~ (.OlOOl + ^Oioioj 

Since the parameter A/" is positive, we have 

Proposition 3.1 (a) Q(°)(Ari°)) = (ACi°\ 0, l^o + aooooAAi°^; A^i^^) is a bifurcation point 
for the approximating system I^^TjEZQ) if Mc? is positive. 

(b) For the double well with well- separation parameter, L, we have that J\fcr\L) > for 
L sufficiently large. 

Proof: We need only check (b). This is easy to see, using the large L approximations of ipo 
and ipi in terms of ipuj, the ground state oi H = —A + V{x), the "single well" operator: 

^Po ~ 2-1/2 (^^(x-L)+V^^(x + L) ) 

i^i ~ 2-1/2 ^ ^^^^ ^^^^ ^ ^) ) . ^3_^^) 

see Proposition 18.11 in section [HI 
Excited state branch 

det{dQj^{Q^'\Af)-Af)) = ^ AfP = (3.12) 

OOIIO + -^O-ooil — Ollll 
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Remark 3.1 For the double well with well-separation parameter, L, we have that J\fi^\L) < 
forL sufficiently large, as can be checked using the approximation h3.11\) . Therefore Q^^^Mi^^) 
is not a bifurcation point of the approximating system liS. 8\) . 

Summary: Assume A/" is sufficiently small. The finite dimensional approximation (13.81) 
predicts a symmetry breaking bifurcation along the nonlinear ground state branch and that 
no bifurcation takes place along the anti-symmetric branch of nonlinear bound states. 



4 Bifurcation / Symmetry breaking analysis of the PDE 

In this section we prove the following 

Theorem 4.1 (Symmetry Breaking for NLS-GP) Consider NLS-GP with hypotheses 
(H2)-(H4). Let akimn be given by l{2.19\) and 

^H'o: i^iyO] = c^oooo " c^iooi ~ 2aioio 

= {i^lgK[i;l]) - i^lgKiij'o]) - 2 i^o^ugK[^oH) > 0- (4.1) 

Assume ^ ^ 

— ^ is sufficiently small. (4.2) 

S[^o,V'i]' 

Then, there exists Mcr > such that 

(i) for any M < Mcr, there is (up to the symmetry u ^ u e^"* ) a unique ground state, Uj^, 
having the same spatial symmetries as the double well. 

(a) M = Mcri "W^™ 0, bifurcation point. For M > Mcr, there are, in a neighborhood 
of M = Mcr, w^'", two branches of solutions: (a) a continuation of the symmetric 
branch, and (b) a new asymmetric branch. 

(Hi) The critical M- value for bifurcation is given approximately by 

^1 i^Q 



H[^o,V'i]' 



Corollary 4.1 Fix a pair of eigenvalues, {Q2j , 'ip2j) , {^2j+i,4'2j+i) of the linear double-well 
potential, Vl{x); see Example \2.1\ For the NLS-GP with double well potential of well- 
separation L, there exists L > 0, such that for all L > L, there is a symmetry breaking 
bifurcation, as described in Theorem\4.1\ with Mcr = Mcr{L] j) ■ 



Remark 4.1 VLi{L) — Qq{L) = 0{e~'^^) forL large. The terms in E[iIjq,iIji]{L) are 0{1). 
Therefore, for the double well potential, Vl{x), the smallness hypothesis of Theorem \4.1\ holds 
provided L is sufficiently large. 
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To prove this theorem we will establish that, under hypotheses fl4.ip - fl4.2l) . the character 
of the solution set (symmetry breaking bifurcation) of the finite dimensional approximation 
f l3.m3.3l) persists for the full (infinite dimensional) problem: 

(fio - r^) Co + aoooolcopco + (aono + aooii) |cipco + aooncjco + {ijjog,TZ{co,ci,r])) = 

(4.3) 

(fii -n)ci + aiiiilcipci + (aioio + Oiooi) |copCi + aioioCoCi + {ipig,n{co, Ci,r])) = 

(4.4) 

(H-n) f] = -P g[F{--co,c,) + 7^(co,Cl,r/)], 7] = Pr] (4.5) 
|co|' + |ci|2 + [ \ri\^ = M. (4.6) 



We analyze this system using the Lyapunov-Schmidt method. The strategy is to solve 
equation (14. 5 p for 77 as a functional of Co,ci and Vt. Then, substituting 77 = rj[cQ,Ci,VL] 
into equations (14. 3p . (14.41) and (14. 6p . we obtain three closed equations, depending on a 
parameter M, for Cq, Ci and f2. This system is a perturbation of the finite dimensional 
(truncated) system: (13.11 13. 2p and (13.30 . We then show that under hypotheses (14. ip - (14. 20 
there is a symmetry breaking bifurcation. Finally, we show that the terms perturbing the 
finite dimensional model have a small and controllable effect on the character of the solution 
set for a range of A/", which includes the bifurcation point. Note that, in the double well 
problem, hypotheses (14. 11) - (14. 2 p are satisfied for L sufficiently large, see Proposition 18. 2[ 

We begin with the following proposition, which characterizes r] = ?7[co, ci, fi]. 

Proposition 4.1 Consider equation for rj. By (H4) we have the following: 

Gap Condition : {ilj - r\ > 2d^ for j = 0, 1 and all r G a{H) \ {^o, ^1} (4.7) 

Then there exists n*, r^, > 0, depending on d^,, such that in the open set 

|co| + |ci| < (4.8) 
||co^/'o + ciV^i + ?7||h2 < n^{d^) 

dist{n,a{H)\{no,ni}) > 4, (4.9) 

the unique solution of (I2.18P is given by the real- analytic mapping: 

{co,Ci,n) t-^ ?7[co,Ci,fi], (4.10) 
defined on the domain given by ^4-^4-^ ■ Moreover there exists > such that: 



II vKcu^ \\H-<a{\co\ + \ci\f (4.11) 

Proof: Consider the map 



A^(0O,01,02) = gK[(pi(p2](p3- 
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By assumptions on the nonlinearity (see section [2]), there exists a constant k > such that 

||iV(0o,01,02)|U2 < k\\MH4MH4h\\H^- (4.12) 

Moreover the map being hnear in each component it is real analytic. 

Let Co, Ci and Q be restricted according the inequalities fl4.8ll4.9p . Equation (12.181) can be 
rewritten in the form 

7] + {H -n)-^pN[cotlJo + ciilJi+7]] = 0. (4.13) 
Since the spectrum of HP is bounded away from Q hj d^, the resolvent: 

is a (complex) analytic map and bounded uniformly, 

\\{H-n)-'P\\L2^H2<p{d;'), (4.14) 

where p{s) oo as s ^ cxd. Consequently the map F : C"^ x {Q E C : dist{Q, a{H) \ 
{^0, ^^i} } > 4} X if2 ^ gjygj^ 

F(co, ci, ^1,7]) = 7] + {H - ^IY'Pn[ co^o + ci^i + ] (4.15) 
is real analytic. Moreover, 

F(0,0,1],0) = 0, ^^^(0, 0, 1], 0) = J. 

Applying the implicit function theorem to equation (I4.13p . we have that there exists n^iVt)^ T*(fi) 
such that whenever |co| + |ci| < r^, and \\CqiI)q + Ciipi + ri\\H2 < equation (14.131) has an 
unique solution: 

which depends analytically on the parameters Cq, Ci, Q. By applying the projection operator 
P to the fl4.13p which commutes with (H—Q)^^ we immediately obtain Pt] = i], i.e. rj G PL^. 

We now show that Ti^j can be chosen independent of f2, satisfying (14. 9p . The implicit 
function theorem can be applied in an open set for which 

D,F(co, ci, fi, r/) = / + (i7 - nY^pDr,N[ coi^o + ci^i + v ] 

is invertible. For this it suffices to have: 

II {H - ny^pDrjNl Co^/Jo + cii)i + ri\ \\h2 < Lip < 1 

A direct application of (I4.12p and (I4.14p shows that 

\\{H - ny'PD^N[ coV'o + ci^i + v] \\h^ < 

3k p{d;') \\cotpo + ci^i + vWl^ (4.16) 
^The trilinearity follows from the implicit bilinearity of K in formulas (|2.16p - (|2.18p . 
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Fix Lip = 3/4. Then, a sufficient condition for invertibility is 

3A; p{d;') IMo + ciV^i + vWm < Lip = 3/4. (4.17) 

which allows us to choose n* = \\J'^~^t^^ independently of VL. 
But, if f HTTIl holds, then, from fH?T6|) . the if^ operator 

(i7-fi)-ipiV[co^o + Ci^i + -] 

is Lipschitz with Lipschitz constant less or equal to Lip = 3/4. The standard contraction 
principle estimate applied to (14.131) gives: 

Mh^ < -^-—m-n)~'PN[c^iJo + CiiJi]\\H^ 
1 — Ltp 

< Md:^) k\Mo + ciMm- (4-18) 
Plugging the above estimate into (14.171) gives: 

llcoV'o + ciipiWn^ + 4:p{d-^) k\\coil)o + ciV'i||^2 < — , 

2VP("* )k 

Since the left hand side is continuous in (co,Ci) G and zero for cq = Ci = one can 
construct > depending only on (i*, k such that the above inequality, hence (14.171) and 
(I4.18p . all hold whenever |co| + |ci| < n. Finally, (14. lip now follows from (|418j).QED 

In particular, for the double well potential we have the following 

Proposition 4.2 Let V = Vl denote the double well potential with well- separation L. There 
exists > 0, e(L*) > and d^:{L^) > such that for L > L^,, we have that for {cq,Ci,Q) 
satisfying dist{Q,a{H) \ {f2o7^i}} } > d^{L^) and \cq\ + \ci\ < e{L^) ?7[co,Ci,f2] is defined 
and analytic and satisfies the hound for some > 0. 

Proof: Since f2o, fii, ipQ and ipi can be controlled, uniformly in L large, via the approxi- 
mations (13. lip , both and r^, in the previous Proposition can be controlled uniformly in L 
large. QED 

Next we study the symmetries of ?7(co, ci, VL) and properties of 7?.(co, ci, rf) which we will 
use in analyzing the equations (I2.16p - (l2.17p . The following result is a direct consequence of 
the symmetries of equation (12.181) and Proposition 14.11 

Proposition 4.3 We have 

Vie'^co, e^^ci, n) = e'^r]{co, Ci, fi), for < ^ < 2%, (4.19) 

r]{co, ci, Q) = r](c^, cI, Q) (4.20) 

in particular 

ri{e''co, ci = 0,n) = e'%{co, ci = 0, fi), (4.21) 

r/(co = 0, e^^ci, n) = e'^(co = 0, Ci, n), (4.22) 
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?7(co,0,f2) is even in xi, ?7(0,Ci,f2) is odd in xi and if Cq, C\ and Q are real valued, then 
r]{cQ,Ci,Q) is real valued. 
In addition 

( ?/^o,^(co,ci,?7) ) = Co fo{ co,ci,n) (4.23) 

( V^i, 7^(co, ci, r])) = h{ Co, ci, ) (4.24) 

where, for any Q <6 <2-n 

fj (e^^co, e^'cu Q ) = fj { co, c,,n), j = 0,l (4.25) 
fj{co,ci,n) = fj{c^,ci,n), J = 0,1 (4.26) 
I/, ( Co,Ci,fi )| < C(|co| + |cl|)^ J = 0,1 (4.27) 

for some constant C > 0. Moreover, both /o and fi can be written as absolutely convergent 
power series: 

f,{co,c,,n)= Yl &L™n(^)cS 4 c^ J = 0,1, (4.28) 

k+l+m+n>A, k—l+m—n=0, m+n=evcn 

where b^ki^ni^) ^''"^ '^^^^ valued when Q is real valued. In particular, if Cq, ci and Vt are real 
valued, then fj {cq,ci,Q ) is real valued and, in polar coordinates, for Cq, Ci 7^ 0, we have 

f,{\co\,\c^\,Ae,^l)= ^Lp(^)e^''^'|coP1cip-, j = 0, 1, (4.29) 

k+m>2, p6Z 

where A6 is the phase difference between ci G C and cq G C 

Proof of Proposition STS We start with (14.191) which clearly implies (14.211) - ( 1422|) . We fix 
Q and suppress dependence on it in subsequent notation. From equation (14.131) we have: 

r/(e*^co,e*'ci) 

= -{H- ny^PN{e''coiJo + e^'ciiJi + r/, e'^Co^o + e'^Ci^i + r/, e'^Co^o + e'^c.^Pi + r]) 
= -{H- ny^Pe'^NicotPo + Ci^Ai + e-'^7], cotpo + c^ipi + e'^^r/, cq^Ao + Ci^Ai + e~'%) 

where we used 

iV(a0i, 6</)2, C03) = a6cAr(0i, 02, h)- (4.30) 

Consequently 

e-'%{e''co, e'^ci) = -{H-Qy^PN[ co4jo+citPi+e~''v, co^Ao+ci^i+e"''r/, co^o+ci^i+e-'^r/] 

which shows that both e~*^?7(e*^co, e*^ci) and ?7(co, ci) satisfy the same equation (14.131) . From 
the uniqueness of the solution proved in Proposition 14.11 we have the relation (14.191) . 

A similar argument (and use of the complex conjugate) leads to ( ]4.20p and to the parities 
of 77(00, 0) and 7^(0, ci). 

To prove (S23D and Km . recall that 

7^(co,Cl,?7(co,Cl,fi)) = iV(co^o + ci^i + 77, co^o + ci^i + ?7, CoV^o + ci^i + 77) 

- N{coipo + cii:i,CQtpo + cii:i,coipo + ci^pi). (4.31) 
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Consider first the case Ci = pi G M. Note that 

( ^ifl','^(co,Pi = 0,?7(co,0)) ) = 0. 

Indeed, for pi = 0, all the functions in the arguments of TZ are even functions (in xi) making 
TZ an even function. Since ipi is odd we get that the above is the integral over the entire 
space of an odd function, i.e. zero. Since ( t/;!, 7^(co, pi, ?7(co, pi)) ) is analytic in pi G M by 
the composition rule, and its Taylor series starts with zero we get (I4.24p for real Ci = pi. To 
extend the result for complex values Ci we use the rotational symmetry of TZ, namely from 
fITO]) . (QUI) and Km we have 

7^ (e^^co, e*^ci, r/(e*^co, e^^Ci, fi)) = e^^7^ (cq, Ci, r/(co, ci, fi)) , < ^ < vr 

hence (14.241) holds for ci = |ci|e^*^ by extending /i via the equality (I4.25p . 

A very similar argument holds for (14.231) . Equation (14.271) follows from the definition of 
7^ and fHTTD . Equation follows from 

We now turn to a proof of the expansions for fj-. (I4.28P and (14.291) . Note first that both 
/o and /i are real analytic in cq, Ci by analyticity of TZ in (I4.23l) - (l4.24p : see (14.3 ip . Note also 
that T] is real analytic by Proposition 14.11 while N is trilinear. Hence, both /q and /i can be 
written in power series of the type (I4.28p . Estimate (I4.27P implies that A; + Z + m + n>4, 
while the rotational invariance (I4.25P implies k — l+m — n = 0. The following parity argument 
shows why m + n hence m — n = I — k and k + / are all even. Assume m + n is odd. Note 
that because of (I4.23p . is the scalar product between an even function (in xi) ipo and 

the term in the power series of TZ in which ipi is repeated m + n times. The latter is an odd 
function (in xi) because ipi is an odd function and it is repeated an odd number of times. 
The scalar product and hence b^i^^ for m + n odd will be zero. A similar argument holds 
for bli,^^, m + n odd. Finally &fc/mn(f^) are real valued when Q is real because they are scalar 
products of real valued functions. 

The form (I4.29P of the power series follows directly from (I4.28P by expressing cq and ci 
in their polar forms: cq = |co|e*^° and Ci = |ci|e*^^, A6 = 6i — 6o, and using that m + n, k + l 
and m — n = —{k — I) are all even. The proof of Proposition 14.31 is now complete. 

4.1 Ground state and excited state branches, pre-bifurcation 

In this section we prove part (i) of Theorem 14.11 as well as a corresponding statement about 
the excited state. In particular, we show that for sufficiently small amplitude, the only 
nonlinear bound state families are those arising via bifurcation from the zero state at the 
eigenvalues Qq and Qi. This is true for general potentials with two bound states. Here, 
however we can determine threshold amplitude, A/'cr, above which the solution set changes. 

A closed system of equations for cq, Ci and Q, parametrized by A/", is obtained upon 
substitution of //[cq, Ci, (Proposition 14. ip into (I4.3ti4.6p . Furthermore, using the structural 
properties (14. 23114.241) of Proposition 14. 3[ we obtain: 

(Qq — Q) Cq + OoOOo|CoPco + (Oono + OoOll) Ici^Cq + O-OOllCiCo + Co/o(Co, Ci, fi) = (4.32) 

(fil - (]) ci + aiiii|cipci + (aioio + aiooi) |copci + aioioCoCi + ci/i(co,ci,fi) = (4.33) 
|coP + |ci|2 +O(|co|' + |ci|2)' = X (4.34) 
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This system of equations is valid for \cq\ + |ci| < r*, independent of L, the distance between 
wells. 

If we choose Ci = 0, then the second equation in the system fl4.32l) is satisfied. In this 
case, a non-trivial solution requires cq 7^ 0. The first equation, (14.331) . after factoring out cq 
becomes 

no-n + ao,oo\co\' + fo{\colO,n) = (4.35) 

where we used (I4.25P to eliminate the phase of the complex quantity cq. Since fl is real (14.351) 
becomes one equation with two real parameters Q, \cq\. Since the right hand side of (14.351) 
vanishes for f2 = f^o and |co| = and since the partial derivative of this function with respect 
to Q, evaluated at this solution, is non-zero, we have by the implicit function theorem that 
there is a unique solution 

n = ng{\co\) = no + aoooolcoP + O(|co|"). (4.36) 

By (14.341) . for small amplitudes, the mapping from |coP + |cip to A/" is invertible. The family 
of solutions 

|co| ^ {\co\e'',\c,\=0,n = ng{\co\)) , eoe[0,27r) 

defined for \cq\ sufficiently small, corresponds to a family of symmetric nonlinear bound 
states, uj^ with ||uAr|li2 = A/", bifurcating from the zero solution at the linear eigenvalue Qq 

UM = ( \co\M^) + r/[|co|,0,fig(|co|](x) ) e''°, 9o G [0,27r) 

n = ngi\co\y, 

see, for example, [20, El]. Since both tpo and ri{\co\,0,Qg) are even (in xi) we infer that uj^f 
is symmetric (even). 

Remark 4.2 A similar result holds for the case cq = leading to the anti- symmetric excited 
state branch. 

Proposition 4.4 For |co| -|- |ci| sufficiently small, these two branches of solutions, are the 
only solutions non-trivial solutions of { 2.1 ). 

Proof: Indeed, suppose the contrary. By local uniqueness of these branches, ensured by the 
implicit function theorem, a solution not already lying on one of these branches must have 
both Co and Ci nonzero. Now, divide the first equation by Cq, the second equation by Ci, and 
subtract the results. Introducing polar coordinates: 

Co = poe*^«, ci = pie''\ AO = 6, - Oq, (4.37) 

we obtain: 

Qi-Qo = aooooPo + («oiio + aooii + aooiie*^^^) pi + fo{ Po,Pi, A6', Q) 

- aiiiiPi - (fliooi + flioio + «ioioe"*^^^) Po - /i ( Po, Pi, ^6^, fi) . (4.38) 

The left hand side is nonzero while the right hand side is continuous, uniformly for Q satis- 
fying (14. 9 p and zero for po = = pi. Equation (14.380 cannot hold for |po| + |pi| < £^ where 
e > is independent of Q. This completes the proof of Proposition 14.41 

Note, however that nothing can prevent (14.380 to hold for larger po and pi possibly leading 
to a third branch of solutions of ( 12.70 . In what follows, we show that this is indeed the case 
and the third branch bifurcates from the ground state one at a critical value of po = p^. 
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4.2 Symmetry breaking bifurcation along the ground state / sym- 
metric branch 

A consequence of the previous section is that there are no bifurcations from the ground state 
branch for sufficiently small amplitude. We now show seek a bifurcating branch of solutions 
to (12.16114.3^ . along which cq ■ ci 7^ 0. As argued just above, along such a new branch one 
must have: 

fio-fi + aooooPo+ («oiio + 00011 + aooiie*^^'')p? + /o(po, Pi, ^6^,^) = (4.39) 
ni-n + annpl+ {aioio + aiooi + aioioe-'^^'^) pI + fi{po,Pi,Ae,n) = (4.40) 

We first derive constraints on A^. Consider the imaginary parts of the two equations 
and use the expansions fl4.29p and the fact that Q is real: 

aooii sm{2Ae)pl + ^ bl^^{n) sm{p2Ae)pf p'r = 

k+m>2, pGZ 

aioio sm{2Ae)pl + Yl ^U(^) Mp^^O)pI' pl"" = 0- 

k+m>2, p& 

Since both left hand sides are convergent series in po, pi, then all their coefficients must be 
zero. Hence sin(2A^) = or, equivalently: 

A^G{o,^,7r,|} (4.41) 

Case 1: G {0,7r}: 

Here, the system f l4.39p -( HI10l) is equivalent with the same system of two real equations with 
three real parameters po > 0, pi > and ^2 : 

Fo(po,Pi,fi) =^ ^0 - + aooooPo + («oiio + 2aooii) P? + /o(po,Pi,fi) = (4.42) 
Fi(po,Pi,fi) =^ f^i - + aiiiiP? + (2aioio + aiooi)Po + /i(Po,Pi,^) = (4.43) 

We shall prove that there is a bifurcation point along the symmetric branch using ( 14.ip -( l4l2|) . 
which depend on discrete eigenvalues and eigenstates of —A + V{x). These properties are 
proved for the double well in section [HI an Appendix on double wells. 

We begin by seeking the point along the ground state branch (pg, 0, ^^(pq)) from which 
a new family of solutions of fl4.42l) - fl4.43l) . parametrized by pi > 0, bifurcates; see (14.361) . 

Recall first that for any po > sufficiently small, Fq (po, 0, f2g(po)) = 0. A candidate for 
a bifurcation point is Pq > for which, in addition, 

Fi(p*,0,fi,(p*)) = (4.44) 

Using (14.11) and (14. 2 p one can check that 

Fi(po, 0, l^g(po)) = fii - fio + (aiooi + 2aioio - aoooo + C'(Po)) Po = (4.45) 
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has a solution: 



Po 



kiooi + 2aioio — flooool 



l + O 



— fig 



kiool + 2aioio — ^ooooP 



(4.46) 



We now show that a new family of solutions bifurcates from the symmetric state at (pq, 0, Qg^p^)). 
This is realized as a unique, one-parameter family of solutions 

Pi ^ iPoiPl),Pl,^asym{pl)) (4.47) 

of the equations: 

Fo(po, Pi, ^) = 0, Fi(po, Pi,^) = (4.48) 

To see this, note that by the preceding discussion we have Fj{pQ, 0, Qgi^p^)) = 0, j = 1, 2. 
Moreover, the Jacobian: 



^^^-^i)(0,p*,l],(p*)) 



2Po(oiooi + 2aioio "~ c^oooo 



d{po,^) 
is nonzero because Po > and 

aiooi + 2aioio - Ooooo + C(po^)) < (4.49) 

since pg solves fl4.45p and f2i — fio > 0. Therefore, by the implicit function theorem, for small 
Pi > 0, there is a unique solution of the system fl4.42l) - fl4.43l) : 

Po = Po(Pi) = Po + ;rT —7^ + C(Po) +C(Pi) (4.50) 

^Po V'^lOOl + ^^1010 — c^oooo / 
^ = ^asymiPl) = ^g{P*o) + pl ( ailll + (2aioiO + QlOOl) '^""" ^"""^^ + C'(pS^) ) + C'(Pl), 

V CLlOOl + ^ctioio — c^oooo / 

(4.51) 

Remark 4.3 (1) Due to equivalence of M and pg + p\ as parameters, for small amplitude, 
we have that symmetry is broken at 

Kr - : (4.52) 

POOOO ~ 0,1001 — ^o-ioiol 

(2) Note also that we have the family of solutions 

e'^ {po{pl)lpO ± Pl^l + ^(Po(Pl), ±Pl, ^^a.ym(pi))) , < < 271, Pl > 0. (4.53) 

Here the ± is present because the phase difference A9 between cq and c\ can be or ir, see 
f l4.4ip and immediately below it. Because po 7^ 7^ pi this branch is neither symmetric nor 
anti- symmetric. Thus, symmetry breaking has taken place. In the case of the double well, 
the ± sign in fl4.53p shows that the bound states on this asymmetric branch tend to localize 
in one of the two wells but not symmetrically in both; see also, f^, fi^ . fl^ . 



19 



Case 2: A9 G {f,f } : 

In both cases the system fl4.39l) - fl4.4Up is equivalent to the same system of two real equations, 
depending on three real parameters po ^ 0, pi > 0, Q : 

Fo(po, Pi, fi) =^ ^lo-^l + flooooPo + aoiioP? + /o(po, Pi, ^) = (4.54) 
Fi(po,pi,f^) =^ 1^1 - + aiiiiPi + aiooiPo + /i(Po,Pi,^) = (4.55) 

As before, in order to have another bifurcation of the symmetric branch it is necessary to 
find a point, ( p^* ,0,Qg{pQ*) ), for which: 

FM\0,ng{p**) = n^-no + {a^ooi-aoooo)Po' + 0{p**^) = 0. (4.56) 

If such a point would exist we will have pg* > pg because Oiooi — aoooo > 2aioio + Oiooi ~ Oqooo 
due to aioio < 0. Hence this bifurcation would occur later along the symmetric branch 
compared to the one obtained in the previous case. Consequently the new branch will be 
unstable because, as we shall see in the next section, it bifurcates from a point where the 
operator already has two negative eigenvalues. 

Moreover, it is often the case (see also the numerical results of section [6]) that the equation 
fl4.56p has no solution due to the wrong sign of the dominant coefficient, i.e. aiooi — Ooooo > 0. 
This can be easily checked, in particular, e.g., for g = —1 and large separation between the 
potential wells, using (13.111) . 



5 Exchange of stability at the bifurcation point 

In this section we consider the dynamic stability of the symmetric and asymmetric waves, 
associated with the branch bifurcating from the zero state at the ground state frequency, 
Qq, of the linear Schrodinger operator —A + V{x)] see figure [TJ The notion of stability with 
which we work is - orbital Lyapunov stability. 

Definition 5.1 The family of nonlinear bound states ^-int . q ^ [0,27r) } is - 

orbitally Lyapunov stable if for every e > there is a S{e) > 0, such that if the initial data 
Uq satisfies 

inf \\uo{-) - *f7(-)e'l//i < S , 

eG[0,27r) 

then for all t ^ 0, the solution u{x,t) satisfies 

inf \\ui;t) - <iln{-)e''\\m < e. 

0e[O,27r) 

In this section we prove the following theorem: 

Theorem 5.1 The symmetric branch is orbitally Lyapunov stable for < po < Pq, or 
equivalently < A/" < J^cr- At the bifurcation point po = Po (A/" = A/'cr), there is a exchange 
of stability from the symmetric branch to the asymmetric branch. In particular, forM > Mcr 
the asymmetric state is stable and the symmetric state is unstable. 
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We summarize basic results on stability and instability. Introduce and L_, real and 
imaginary parts, respectively, of the linearized operators about "^q: 



L 



'+ 



iH-n)-+D^N[^n]{-) 



L 



-1 



(5.1) 



By (122D and ([23D, L^^a = 0. 



We state a special case of known results on stability and instability, directly applicable 
to the symmetric branch which bifurcates from the zero state at the ground state frequency 



(1) (Stability) Suppose L+ has exactly one negative eigenvalue and L_ is non-negative. 
Assume that 



Then, is orbitally stable. 

(2) (Instability) Suppose L_ is non-negative. > 2 then the linearized dynam- 

ics about \l/f7 has spatially localized solution which is exponentially growing in time. 
Moreover, is not orbitally stable. 

First we claim that along the branch of symmetric solutions, bifurcating from the zero 
solution at frequency f2o; the hypothesis on L„ holds. To see that the operator 
is always non-negative, consider L„[\E'f7(,] = L_[0] = —A + V — VIq. Clearly, L^[Q] is a 
non-negative operator because VLq is the lowest eigenvalue of —A -|- V . Since clearly we have 
L.^^Q = 0, G spec{L_\^Q\). Since the lowest eigenvalue is necessarily simple, by continuity 
there cannot be any negative eigenvalues for VL sufficiently close to VLq. Finally, if for some 
fi, L_ has a negative eigevalue, then by continuity there would be an for which L_\^q^ 
would have a double eigenvalue at zero and no negative spectrum. But this contradicts that 
the ground state is simple. Therefore, it is the quantity n_(L+), which controls whether or 
not \E'n is stable. 

Next we discuss the slope condition (15.21) . It is clear from the construction of the 
branch Vt i— > \E'f7 that (15.21) holds for Vt near VLq. Suppose now that J|\EfQp = 0. 
Then, (\E'f7, d^^^o) = 0. As shown below, has exactly one negative eigenvalue for 
fl sufficiently near Qq. It follows that > on {\E'n}"'" [25l [26] . Therefore, we have 



iLldn-^n,Lldn-^n) = {L_,dn'^n, dn^n) = (^Q,9f,^n) = 0. Therefore, Lldn^n = 0, im- 



plying = L_|_9f7\E'n = 0, which is a contradiction. It follows that (15. 2p holds so long as 
L_|_ > on {\E'n}"'" and when (15.21) first fails, it does so due to a non-trivial element of the 
nullspace of L+. 

Therefore is stable so long as n_{L^) does not increase. We shall now show that for 
Af < Men '^-{L+\^q\) = 1 but that along the symmetric branch for Af > Ncr ^^-(-^+[^^1) = 
2. Furthermore, we show that along the bifurcating asymmetric branch, the hypotheses of 
Theorem 15.21 ensuring stability hold. 



of -A + V. 



Theorem 5.2 JM\M\^ 




(5.2) 
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Remark 5.1 For simplicity we have considered the most important case, where there is 
a transition from dynamical stability to dynamical instability along the symmetric branch, 
bifurcating from the ground state of H . However, our analysis which actually shows that 
along any symmetric branch, associated with any of the eigenvalues, ^2j,j > of H , there 

is a critical N" = Mcrii), such that as M is increased through Mcr{j), the number of 

negative eigenvalues of the linearization about the symmetric state along the j*^ symmetric 
branch increases by one. By the results in JTTl 0, [7^ . this has implications for the number 
of unstable modes of higher order (j > 1) symmetric states. 

Consider the spectral problem for L+ = L4.[\E'q]: 

L+[^n]0 = /i0 (5.3) 

We now formulate a Lyapunov-Schmidt reduction of f l5.3p and then relate it to our formu- 
lation for nonlinear bound states. We first decompose relative to the states ipQ, ipi and 
their orthogonal complement: 

(f) = aoiJo + aiipi + ^, (^/;j-,^) = 0, J = 0,1 

Projecting (15. Sp onto ipo, ipi and onto the range of P we obtain the system: 

{ijjo,L+['^n\{aoipo + aiipi + ^) ) = jjao (5.4) 

{^i,L+[^Q]{ao^po + ai^i + ) = /lai (5.5) 

(ff-fi)e + A.iV[^n](«oV^o + «i^i + = (5.6) 
The last equation can be rewritten in the form: 

J + (/7-fi-^)-ipZ}„iV[*f,]]e = -{H~n-^Ji)~^PD^N[^>n]{a,^^ + aM (5.7) 

The operator on the right hand side of (15.71) is essentially the Jacobian studied in the proof of 
Proposition 14. 11 evaluated at 17 + /i. Hence, by the proof of Proposition 14.11 if f2 + /i satisfies 
(USD and ll^t^llna < AC, then the operator I+{H -Vt- ii)"^PDuN[^!q\ is invertible on if^ 
and (15.71) has a unique solution 

= g[/i, ^'n](ao^o + ai^i) (5.8) 
= -{I + [H-n- fir^PD^N[m^])-\H -n- fi)-'PD^N[<i/n]Mo + «iV^i) 
= C [(IpoI + |pi|)^] [ ao^Ao + ai^i] . 

The last relation follows from DuN[ip] being a quadratic form in = poV^o+Pi'^i + C^((|po| + 
IPil)')- 

Substitution of the expression for ,^ as a functional of aj into (15.41) and (15.51) we get a 
closed system of two real equations: 

(fio-fi)ao + {^o,DuN[^n] (/ + g[p,^n]) (aoV^o + ai^Ai)) = P 

(fii-fi)ai + (7/;i,D„A^[^'f^] (/ + g[/i,^'n]) (ao^o + aiV^i)) = /i ai (5.9) 

The system (15.91) is the Lyapunov Schmidt reduction of the linear eigenvalue problem for 
with eigenvalue parameter fi. Our next step will be to write it in a form, relating it to the 
linearization of the Lyapunov Schmidt reduction of the nonlinear problem. 
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Remark 5.2 For W'^nWn^ < f^*, the above system is equivalent to the eigenvalue problem 
for the operator L^[^q] with eigenvalue parameter /i as long as \4-9{ ) holds with Q replaced by 
fl + fi. This restriction on the spectral parameter, fi, is in fact very mild and has no impact 
on the analysis. This is because we are primarily interested in fi near zero, as we are are 
interested in detecting the crossing of an eigenvalue of from positive to negative reals as 
M is increased beyond some Mcr- Values of fi for which Iji4-^ does not hold, do not play a 
role in any change of index, n_{L^). 

First rewrite (15. 9p as 

{no-n- ^)ao + {^0, D^N[^n] (/ + g[0, ^n]) (aoV^o + ai^i)) (5.10) 
+ {iJo,D^N[^n] Ag[^,^n] (aoV^o + ttiV^i)) = 

+ {^,,DuN[<i'n] {I + Q[0,<i!n\) (ao^o + "i^^i)) 

+ (V'l, DuN[^n] AQ[fi, ^n] {ao^o + ai^Ai)) = 0. (5.11) 

Here, 

AQ[fi,<iIn] = Qk^d - Q[0,*d- (5.12) 
Note that terms involving AQ in (15. 10115. lip are of size 0[{pq + pDfiaj]. 

Proposition 5.1 

Q[0,'^n]{ao'(po + ai^i) = dp^r][po, pi,n] ao + dp^r][po, pi,n] ai (5.13) 
Proof of Proposition 15. il Recall that rj satisfies 

F{po,pi,n,ri) = r] + {H-ny'pN[poiJo + Pi^i+v] = 0, (5.14) 
Differentiation with respect to pj, j = 0, 1 yields 

(j + {H-ny'PD^N[^n]) dp^T] = -{H-n)-'PD^N[^n]^Pj, (5.15) 

where 

*Q = poiPo + pi'4'i + 'n[po,pi,^- 

Thus, 

ri = Q[0,^n] i^j, (5.16) 

from which Proposition 15.11 follows. 

We now use Proposition 15. II to rewrite the first inner products in equations fl5.10p - fl5.lip . 

For A; = 0, 1 



^{ipk, D^N[pQipQ + piipi + ri\{ipj + dp^vi))i 

j=0 

1 ^ 
1 

^ dp^ {ipk, N[poiJo + pi^pi]) aj + dp^ [ pkfkipo, Pu^) ] , (5.17) 

j=0 
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where iV[^n] = A^[po^o + Pi^i] + ^; see equations 02. 16112. 18p . fl4.23ll4.24l) . Therefore, the 
Lyapunov-Schmidt reduction of the eigenvalue problem for becomes 

(fio-^-At)tto + E]=o ^Pj (V^o,A^[Po^o + PiV^i]) ttj + dpj [ Po/o(Po,Pi,fi) ] (5.18) 
+ {iP^,DuN[^^] Ag[/i,^t^] (ao^Ao + aiV^i)) = 

(f]i-fi-/i)ai + Y.]=o {ipi, N[poipo + piipi] ) aj + dp^ [ pi/i(po, Pi, ^^) ] 

+ (tAi, Z^„iV[*n] Ag[/i, («o^o + tti^i)) = 0. (5.19) 



M - /i + C(/i) ] ( !° 1 = f !! ) , (5.20) 



This can be written succinctly in matrix form as 

ai j y 

where 

M = M[17,po,Pi] 

fio - + SaooooPo + («oiio + 2aooii)p? + dpg{pofo) 2(aoiio + 2aooii)poPi + dp^{pofo) 

2(2aioio + aiooi)poPi + i9po(Pi/i) (^i - ^) + Saiiiipj + (2aioio + aiooi)pl + dp,{pifi) 

(5.21) 



and 

Note that 



C{fi)im = {^u D^N[^^]^Q[ii, /, m = 0, 1. (5.22) 



C(/i = 0) = 0. (5.23) 

Recall that /i is the spectral parameter for the eigenvalue problem L+, fl5.3p and we are 
interested in n_(L+[^t7]), the number of negative eigenvalues along a family of nonlinear 
bound states ^ By Theorem 15.21 n_ (L_|_) determines the stability or instability of a 
particular state. This question has now been mapped to the problem of following the roots 
of 

D{p, po, Pi) = det(/i/ - M - C(p)) = 0, (5.24) 

where po and pi are parameters along the different branches of nonlinear bounds states. 
Since C(p), defined in fl5.22p is small for small amplitude nonlinear bound states, we expect 
the roots, p, to be small perturbations of the eigenvalues of the matrix M. We study 
these roots along the symmetric (M = M(f2g(po), po, 0)) and asymmetric branch (M = 
M {VLasymipi) , Pq{,Pi) , Pi)) usiug the implicit function theorem. 

Symmetric branch: 

Along the symmetric branch: 

Pi = 0, Po > 0, = fig = r^o + aooooPo + C'(Po)' = PoV^o + '7(po, 0, fi) = symmetric. 

Thus, D = D{ji,po). Moreover, the system fl5.20p is diagonal. This is because Q, and hence 
AQ, preserve parity at a symmetric see their definitions. Therefore Coi = = Cio, 
each the scalar product of an even and an odd function. Moreover from fl4.29p we get: 
|^(po,0,fi) =0, J =0,1. 
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Therefore, the matrix fil — M — C(/i) is diagonal and is an eigenvalue of -^^+['?/'Cg(po)] if 
and only if is a root of either 

^o(/W, Po) = At - ^oo(Po) - Coo(p, Po) = (5.25) 

or 

Pi(/i,Po) = /i-Mn(po) -Cii(p,Po) =0 (5.26) 
Both Pq and Pi are analytic in fi and po and it is easy to check that 

Po(0,0) =0, 9^Po(0,0) = l 

and 

Pi (1^1 - 1^0, 0) = 0, d^Piin^ - no, 0) = 1. 

Formally differentiating fl5.25p or (15.261) with respect to po gives: 

''p'^' ~ 1 - d,c,, • ^^-^^^ 

By the implicit function theorem (15.251) and (I5.26P define, respectively, po and pi as smooth 
functions of p provided 

\d^C,,\<l, J = 0,1 (5.28) 

A direct calculation using (15. 8p and estimates (I4.12p . (14.140 shows that in the regime of 
interest: VL satisfying (14.90 . it suffices to have 

W'^q.Wh^ < n, (9max(||^o||H2, UiWh^))'^ (5-29) 

where n* is given by Proposition 14.11 The latter can be reduced to an estimate on po via the 
above definition of and (14.180 as in the end of the proof of Proposition 14. 1[ 

Therefore, under conditions (14. 9 p and (15.290 . we have a unique solution po, respectively 
Pi, of (15.250 . respectively (15.260 . Moreover, the two solutions are analytic in po and, for 
small Po, we have the following estimates: 

Po = 2aooooPo + (^(Po) < (5.30) 
pi = fil - 1^0 + C(pg) > 0, (5.31) 

where we used aoooo = g{4'oy ^bPol) < 0' Pi(Po = 0) = f2i — f2o > 0. 

We claim that pi changes sign for the first time at po = p^. Indeed, by continuity, the 
sign can only change when pi = 0, i.e. when (I5.26P has a solution of the form (0,po). Since 
Cii(0,po) = 0, see ([S2SD, fl5^ becomes 

= Mii(po) = fil - fig(po) + (2aioio + aiooi)Po + /i(Po, 0, fig) = Pi(po, 0, fi9(po)); 

see (I5.2ip and note that pi = 0. But this equation is the same as (14.440 . which determines Po, 
then bifurcation point. Thus, as expected, -D(p, po) = has a root pi(Po) = or equivalently 
L+ has a zero eigenvalue at the bifurcation point. Note that the associated null eigenfunction 
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has odd parity in one space dimension, and is more generally, non-symmetric and changes 
sign. 

To see that fii{po) changes sign at po = Po differentiate fl5.26p with respect to po at 
Po = Po a^cl obtain from (15.271) that 



^""^^ - 1 - d,Cu ^ 



This follows because the denominator is positive, by fl5.28p . while direct calculation gives for 
the numerator: 

dpoMiiipl) + dp^Cuipl) = 2pS (oiooi + 2aioio - aoooo + C'(pS^)) < 

see (14.491) . Therefore pi becomes negative for po > Pg at least when |po — Po| is small enough. 

In conclusion, L^[Qg{po)] has exactly one negative eigenvalue for < po < Po and two 
negative eigenvalues for po > Pg and |po — Po| small. Therefore, following the criteria of 
[25| [26l [71 [m [6l [12] , the symmetric branch is stable for < po < Pg and becomes unstable 
past the bifurcation point. 

Asymmetric branch: Stability for > Mcr 

Finally, we study the behavior of eigenvalue problem (15.201) on the asymmetric branch: 

< pi < 1 

/ \ * , Pi / '^0110 + 2aooil — 0,1111 , *2^^ , /o^ 4n /c- Qo^ 

Po = Po(Pi) = Po + TTT -7^ + ^(Po ) + C'(Pi) (5.32) 

zPq \aiooi + ^Oioio ~ ooooo / 

on t \ n ^ *^ , 2 , /o , x Ooiio + 2aooii - aim , *2^^ 

il = ilasymiPl) = i^g(Po) + Pi O-UU + (2aioiO + aiooi) — h 0{Pq ) 

\ OlOOl + ^OlOlO ~ ooooo / 

+0{pt), (5.33) 

= Po(Pl)^0 + PlV^l +'7(Po(Pl),Pl,^as3/m(Pl)) 

The eigenvalues will be given by the zeros of the real valued function 

D{p, pi) = det(p/ - M(pi) - C(p, pi)), (5.34) 

which is analytic in p and pi for Q + p satisfying (14. 9 p and pi small. Note that at pi = 
we are still on the symmetric branch at the bifurcation point Po = Po- Hence, the matrix is 
diagonal and 

D(p,0)=Po(p,pS)A(p,pS), (5.35) 
where Pj, j = 0,1 are defined in (I5.25p - (l5.26p . In the previous subsection we showed that 
each Pj(-, Pq) has exactly one zero, pj, on the interval — oo < p < d^, — flgipo) > 0. The zeros 
were simple, by our implicit function theorem application in which, 

dpP,ip„ pS) = 1 - dpC,, > 0, (5.36) 

see (15.281) . In addition one can easily deduce that lim^^_oo -P; (/i, Po) = — oo by using the 
definitions (15.221) . (15.121) and the fact that \\{H — f2 — p)^^\\l'^^h^ ^ — which implies 
\\Q[p,^n]\\H^^H^'"^-^'^ 0. 
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Consequently D{-,0) has exactly two simple zeros /iq < and /ii = on the interval 
— oo < < {—d^: — VLg{pl))/2 > 0, which are both simple and lim^_,_oo -D(/i, 0) = oo. It is 
well known, and a consequence of continuity arguments and of the implicit function theorem, 
that the previous statement is stable with respect to small perturbations. More precisely, 
there exists £ > such that whenever |pi| < e, D{-, pi) has exactly two zeros po{pi) < and 
/ii(pi) on the interval —oo < p < (— (i* — fig(po))/2 > 0, which are both simple and analytic 
in pi. 

Since we are interested in n_(L_|_), the number of negative eigenvalues of L_|_, we still 
need to determine the sign of pi{pi). In what follows we will show that its derivatives satisfy 

5p,/ii(0)=0, 9>i(0)>0. (5.37) 

We can then conclude that for < pi ^ 1, pi(pi) > 0, and has exactly one (simple) 
negative eigenvalue, po(Pi)- Therefore, the asymmetric branch is stable. 
We now prove (15.371) . By differentiating 

D(pi(pi),pi) = (5.38) 

once with respect to pi at pi = we get 

a^D(0,0)9,,pi(0) + 9p,D(0,0) = 

Using (I5.35P we obtain 

a^D(0,0) = Po(0,Po)5^Pi(/ii = 0,p*) > (5.39) 
where we used (lOHl) and that Po(0, Po) = -^oo(pS) > 0- Using (Oil) and we obtain 

9p,L>(0, 0) = ^^^^^{pi = 0) = det 10 + det 01, where (5.40) 
opi 



det ij = the determinant evaluated at pi = of the matrix obtained from M by differen- 
tiating the first row i times, respectively the second row j times, det ij can be evaluated 
using (gHHD, (ICTD . and fl533D . 

Note that the second row of det 10 is zero and therefore det 10 = 0. Furthermore, det 01 is 
zero because its second column is zero. Therefore, by (I5.40p we have (9p^pi(0) = 0.. 

We now calculate 9^^pi(pi = 0). Differentiate fl5.38p twice with respect to pi at pi = 
and use 9p^pi(0) = to obtain: 

d,D{0,0)dlpi{0) + dlD{0,0)=0. 

which implies, by (15.391) 

sign(9>i(0)) = -sign(92^Z}(0,0)). 
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But, as before, (15.341) and fl5.23p imply 



D(0,0) = ^ : M O) = det20 + 2detll + det02 < 0. 

The last inequality is a consequence of the following argument. First, det20 = 0, since its 
second row zero. A direct calculation using the definition of M and relations fl5.32p show: 

detll = -4(aoiio + 2aooii)(2aioio + oiooOpS^ + C(pS^) 
det02 = SaooooaiiiiPo' + (^(pS") 

Note that in the limit of large well-separation hmit (L >> 1), all coefficients akimn = o,kimn{.L) 
converge to the same value ga^ < 0. This implies 

2 det 11 + det 02 = {-Ug'^a^ + 0{e^^^))pf + C(p*^) < 0. 

Therefore, d'l /ii(0) > and the proof of Theorem 15. II is now complete. 



6 Numerical study of symmetry breaking 

Symmetry breaking bifurcation for fixed well-separation, L 

In this section we numerically compute the bifurcation diagram for the lowest energy 
nonlinear bound state branch for NLS-GP (12. ip and compare these results to the predictions 
of the finite dimensional approximation Eqs. (13.81 Specifically, we numerically compute the 
bifurcation structure of Eq. (I2.ip for a double- well potential, Vl{x), of the form: 

V{x) = V^ ^=exp — ^ +^=exp • 6.1 

The potential for Vq = — 1 , s = 1 and L = 6 has two discrete eigenvalues f2o = —0.1616 
and = —0.12 and a continuous spectral part for > 0. The linear eigenstates can 
also be obtained and used to numerically compute the coefficients of the finite dimensional 
decomposition of Eqs. (13. 8p as aoooo = —0.09397, aim = —0.10375, aooii = aioio = c^iooi = 
ooiio = —0.08836 (for g = —1). Then, using (I3.10p . we can compute the approximate 
threshold in A/" for bifurcation of an asymmetric branch (and the destabilization of the 
symmetric one): 

Kr ~ Af^? = 0.24331, Qcr ~ = ^0 + aoooo Af^? = -0.18447. 

We expect good agreement because the values of s and L suggest the regime of large L, 
where our rigorous theory holds. 

Using numerical fixed-point iterations (in particular Newton's method), we have obtain 
the branches of the nonlinear eigenvalue problem (12. 7p . To study the stability of a solution, 
uq, of (12.70 . consider the evolution of a small perturbation of it: 

(6.2) 



u 



-int 



Uo(X 



p{x] 



q{x) 



28 



Keeping only linear terms in p, g, we obtain a linear evolution equation, whose normal 
modes satisfy a linear eigenvalue problem with spectral parameter, which we denote by A 
and eigenvector {p{x) , q{x))'^ . 

Our computations for the simplest case of the cubic nonlinearity with Klipip] = i/jip 
are shown in Figure [3] (for g{x) = —1). In particular, the top subplot of panel (a) shows 
the full numerical results by thin lines (solid for the symmetric solution, dashed for the 
bifurcating asymmetric and dash-dotted for the anti-symmetric one) and compares them with 
the predictions based on the finite dimensional truncation, (13.81) shown by the corresponding 
thick lines. The approximate threshold values J\fcr and flcr are found numerically to be 
fii? ^ -0.1835, Af^r' ~ 0.229. This suggests a relative error in its evaluation by the finite- 
dimensional reduction of less than 1%. This critical point is indicated by a solid vertical 
(black) line in panel (a). For > f2c?, there exist two branches in the problem, namely 
the one that bifurcates from the symmetric linear state (this branch exists for < Qq) and 
the one that bifurcates from the anti-symmetric linear state (and, hence, exists for Q < Qi). 
For f2 < f2ir\ the symmetric branch becomes unstable due to a real eigenvalue (see bottom 
subplot of panel (a)), signalling the emergence of a new branch, namely the asymmetric one. 
All three branches are shown for Q = —0.25 (indicated by dashed vertical (black) line in 
panel (a)) in panel (b) and their corresponding linearization spectrum (Ar, Aj) is shown for 
the eigenvalues \ = Xr + iXi. 

Symmetry breaking threshold, McriL) as L varies 

We now investigate the limits of validity of N'cr\L) as an approximation to AfcriL) by 
varying the distance L between the potential wells (16.11) . For L large, J\fcr \ given by equation 
(I3.10p . is close to the actual Mcr{L), the exact threshold. In this case the eigenvalues of 
—(9^ -|- Vl[x), 0.q{L) and VLi{L), are close to each other; see Remark HI Therefore, the 
bifurcation occurs for small M and one is in the regime of validity of Theorem 14. 1[ In figure 
m we display a comparison between the estimate for Mcr based on the finite dimensional 
truncation, Ncr\ and the actual Mcr- For large L the two values are close to each other. 
As L is decreased the wells approach one another and eventually, at L = 0, merge to form 
a single well potential. As L is decreased, the eigenvalues of the linear bound states VLq{L) 
and fii(L) move farther apart. For some value of L, L^, the eigenvalue of the excited state, 
r2i(L), merges at = 0, into the continuous spectrum. For L < Ld the estimate Mc? is 
not correct. In fact, J\fcr\L) —>■ oo, while the actual value of J\fcr{L) appears to be remain 
finite. In Figure we observe that for L < 2, A/"^"-* and A/'cr diverge from one another 
and eventually the approximation Mcr\L) tends to infinity, while the actual Mcr{L) remains 
finite. Moreover, in FigureHb we show a bifurcation diagram for small L in which the discrete 
(excited state) eigenvalue of —d1 + V/^, fii, does not exist, and yet there exists a symmetry 
breaking point Mcr- 

More general nonlinearities 

To simplify the analysis, we assumed a cubic nonlinearity in NLS-GP. The analogue of 
the finite-dimensional approximation (13.81) can be derived, for more general nonlinearities, 
by the same method. In this section we present numerical computations for general power 
law nonlinearities such as Ky^pil)] = {ipipY and observe similar phenomena to the cubic case 
p = 1. This is illustrated e.g. in Figure [5l presenting our numerical results for the quintic 
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Figure 3: (Color Online) The figure shows the typical numerical bifurcation results for the 
cubic case and their comparison with the finite dimensional analysis of Section [3l Panel 
(a) shows the bifurcation diagram in the top subplot and the relevant real eigenvalues in 
the bottom subplot. In the top, the solid (blue) line represents the symmetric branch, the 
dash-dotted (green) line the anti-symmetric branch, while the dashed (red) line represents 
the bifurcating asymmetric branch. The thin lines indicate the numerical findings, while the 
thick ones show the corresponding finite-dimensional, weakly nonlinear predictions. The solid 
vertical (black) line indicates the critical point (of Q ~ —0.1835) obtained numerically. The 
dashed vertical (black) line is a guide to the eye for the case with Q = —0.25, whose detailed 
results are shown in panel (b). The bottom subplot of panel (a) shows the real eigenvalue (as 
a function of fl) of the symmetric branch that becomes unstable for Q < —0.1835. Panel (b) 
shows using the same symbolism as panel (a) the symmetric (left), anti-symmetric (middle) 
and asymmetric (right) branches and their linearization eigenvalues (bottom subplots) for 
Q = —0.25. The potential is shown by a dotted black line. 



case of p = 2 (the relevant curves are analogous to those of Figure [3]) • It can be observed that 
the higher order case possesses a similar bifurcation diagram as the cubic case. However, the 
critical point for the emergence of the asymmetric branch is now shifted to n'S^ ^ -0.1725, 
i.e., considerably closer to the linear limit. In fact, we have also examined the septic case of 
p = 3, finding that the relevant critical point is further shifted in the latter to fie? = —0.168. 
This can be easily understood cLS CcLSGS with higher p are well-known to be more prone to 
collapse-type instabilities (see e.g. [23]). It may be an interesting separate venture to identify 
fie? as a function of p, and possibly obtain a such that Vfi < fio, the symmetric branch is 
unstable. We also note in passing that bifurcation diagrams for higher values of p may also 
bear additional (to the shift in fii?) differences from the cubic case; one such example in 
Figure [5] is given by the presence of a linear instability (due to a complex eigenvalue quartet 
emerging for fi < —0.224) for the anti-symmetric branch. The latter was found to be linearly 
stable in the cubic case of Fig. [3l 

Nonlocal nonlinearities 
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Figure 4: (Color Online) The figure demonstrates the validity oiNcr {L) as an approximation 
to McriL). Panel (a) compares the linear finite dimensional estimation for the bifurcation 
point Mcr\L) and the actual numerical bifurcation point Mcr- The computations are for the 
double well potential (16.11) Vq = — 1 and s = \ and cubic nonlinearity. The curve Afcr{L) 

is marked by a solid (black) line and the curve Mcr\L) is marked by a dotted (blue) line. 
Panel (b) shows a numerical bifurcation diagram for the double well potential (16.11) Vq = — 1, 
s = \ and L = 1.3. The bifurcation point Afcr is marked by a (red) circle. For Af < Mcr the 
ground state marked by a thick (blue) solid line is stable. For M > Mcr the ground state 
is unstable and marked by a thick (blue) dashed line. The stable asymmetric state which 
appears for M > Mcr is marked by a thin (red) solid line. The unstable antisymmetric state 
(f2^) is marked by a thin (light green) dashed line. The point M for which the antisymmetric 
state appears in the discrete spectrum is marked by a (black) square. Notice that in this 
bifurcation diagram there is also a bifurcation from the antisymmetric branch. The state 
which bifurcates from the antisymmetric state is marked by a (dark green) thin dotted line. 



Finally, we consider the case of nonlocal nonlinearities, depending on a parameter e, the 
range of the nonlocal interaction. In particular, consider the case of a non-local nonlinearity 
of the form: 

/oo 
IC{x - y),p{y)^{y)dy, (6.3) 
-oo 

where 

lC{x-y) = ^y^. (6.4) 

Here, e > is a parameter controlling the range of the non-local interaction. As e tends to 
0, K,{x — y) ^ 5{x — y) and we recover the "local" cubic limit, limit. The form of the finite 
dimensional reduction, (13. 8p . does not change; the only modification is that the coefficients 
O'kimn are now functions of the range of the interaction e. The dependence of the coefficients, 
O'kimn on e is displayed in panel (a) of Fig. [61 The solid (blue) line shows laooool; the dashed 
(green) one corresponds to |aiiii|, the dash-dotted (red) one to |aiooi| = looiiol (due to 
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Figure 5: Same as Figure [3] but for the quintic nonlinearity. This serves to illustrate the 
analogies between the bifurcation pictures but also their differences (shifted critical point 
and also partial instability of the anti-symmetric branch). 

symmetry), while the thick solid (black) one to |aoioi| = kooiil = Ic^ioiol- Notice in the 
inset how the coefficients asymptote smoothly to their "local" limit. Additionally, note the 
expected asymptotic relation aiooi = ctooii- Also note the significant (decaying) dependence 
of the relevant coefficients on the range of the interaction. The nature of this dependence 
indicates that while the character of the bifurcation may be the same as in the case of local 
nonlinearities, its details (such as the location of the critical points) depend sensitively on 
the range of the non-local interaction. This is illustrated in panel (b) for the specific case 
of e = 5. In this panel (which is analogous to panel (a) of Figure [3], but for the non-local 
case) the critical point for emergence of the asymmetric branch/instability of the symmetric 
branch is shifted to fie? = —0.2466 (and the corresponding A^r = 1.4353) in comparison to 
the numerically obtained value of fii? ~ —0.256; the relative error in the identification of 
the critical point (by the finite-dimensional reduction) is in this case of the order of 3.7%, 
which can be attributed to the more strongly nonlinear (i.e., occurring for higher value of 
nature of the bifurcation. However, as the finite-dimensional approximation still yields 
a reliable estimate for the location of the critical point, in panel (c) we use it to obtain an 
approximation to the location of the critical point (fii? , AClr ^ ) as a function of the non-locality 
parameter e. 




7 Concluding remarks 

We have studied the spontaneous symmetry breaking for a large class of NLS-GP equations, 
with double-well potentials. Our analysis of the symmetry breaking bifurcation and the 
exchange of stability is based on an expansion, which to leading order in amplitude, is a 
superposition of a symmetric - antisymmetric pair of eigenstates of the linear Hamiltonian, 
H, whose energies are separated (gap condition ( 14. 7p ) from all other spectra of H. This 
gap condition holds for sufficiently large L but breaks down as L decreases. Nevertheless, 
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numerical studies show the existence of a finite threshold for symmetry breaking. A theory 
encompassing this phenomenon is of interest and is currently under investigation. 

8 Appendix - Double wells 

In this discussion, we are going to follow the analysis of |B]. Consider a (single well) real 
valued potential Vo{x) on M" such that Vq{x) & U + for all 1 < r < g where q > 
max(n/2, 2) for n 7^ 4, g > 2 for n = 4. Then, multiplication by Vq is a compact operator 
from to and 

Ho = -A + voix) 

is a self adjoint operator on with domain H^. 
Consider now the double well potential: 

Vl = TlVoT^l + RTlVqT.lR 
where Tl and R are the unitary operators: 

TlQ {Xi,X2, ... ,Xn) = g{Xi + L,X2, . . . ,Xn) 

Rg {Xi,X2, ... ,Xn) = g{-xi,X2, . . . ,Xn) 
and the self adjoint operator: 

Hl = -A + Vl{x) 

Proposition 8.1 Assume that u < is a nondegenerate eigenvalue of Hq separated from 
the rest of the spectrum of Hq by a distance greater than 2(i*. Denote by ip^^ its corresponding 
e-vector, \\'4'u)\\l'^ = 1- Then there exists Lq > such that for L > Lq the following are true: 

(i) Hl has exactly two eigenvalues ^lo{L) and fii(L) nearer to uj than 2d^. Moreover 
limi^oo ^j{.L) =uj, j = 0, 1. 

(a) One can choose the normalized eigenvectors ipj{L), \\iljj{L)\\L2 = 1, corresponding to 
the e-values Qj{L), j = 0, 1 such that they satisfy: 

lim \\ijj{L) - {TLiJu. + {-iyRTLi'^)/V2\\H2 = 0, j = 0, 1. 

(Hi) If Pj" are the orthogonal projections in onto ipj{L), j = 0, 1 and Pl = Id — Pq — P^ 
then there exists d > independent of L such that: 

\\{IIl — Q)^^ Pl\\l2^jj2 > (i, for all L > Lq and — uj\\ < d^. 

Proof: For (i) we refer the reader to [S] . The convergence in (ii) has also been proved 
there. The H"^ convergence follows from the following compactness argument. Let: 

^} = nLij,{L), j = 0,l 
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where ul is such that ||'?/'j'||H2 = 1, j = 0,1. From the eigenvector equations: {H^ — 
Q{L))iIj^ = 0, where we dropped the index j = 0, 1 and the convergence Q{L) — > u, see part 
(i), we get 

hm II (-A -u + VL)tp^\\L2 = 0. (8.1) 

L— >oo 

Denote: 

gL = i-A-co)ij'^ eL\ (8.2) 

Since —A — u : t— is bounded there exists a constant C > independent of L such 
that 

||5'l||l2 < c. 

Since u; < 0, — A — u; : has a continuous inverse then (18. ip is equivalent to: 

(7L + yL(-A-^)-i(?i^0, mL\ 

By expanding Vl we get 

9l + Tii;o(-A - c^)-^T_i^?L + RTlVo{-A - uoY^T^lRQl ^ 0. (8.3) 

But vq{—A — uj)^^ : ^ L'^ is compact while the translation and reflection operators are 
unitary. These and the uniform boundedness of ql lead to the existence oi ijj E L"^ and 
ip E L"^ and a subsequence oi Ql, which we will redenote hy ql, such that 

lim ||t;o(-A - uj)-'T^L9L - ^lU^ = and lim ||t;o(-A - uY^T^LRgi - ^Wl^ = 0. (8.4) 

L—*oo L— >oo 

By plugging in (18. 3p and multiplying to the left by T^l we get 

lim \\T^L9L + i^ + RT2Li'\\L^ =0. 
L— >oo 

But RT2L'ip converges weakly to zero, hence T^lOl converges weakly to —tp. By plugging 
now in (18.41) and using compactness we get: 

i/j + voi-A - cj)"V = 0. 

The latter shows that (—A — uj)~^tlj is an eigenvector of —A + vq corresponding to the 
eigenvalue u. By nondegeneracy of uj we get 

iP = -n{-A-uj)tlj^, (8.5) 

where n is a constant. A similar argument shows 

^P = -n{-A-u;)^P^, (8.6) 

where n is a constant. 

Combining fl8.1l) - fl8.6p we get 

lim ||(-A - a;)(^/>^ - uTl^P^ - ni?Ti^^)|U2 = (8.7) 

L—>oo 
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which by the continuity of (—A — u) ^ : h-> imphes 

hm W^p^ - riTii^^ - hRTLtpJln^ = 

Using now that ||^/^^||j:/2 = 1 and that the rescaled ipj such that it has norm 1 in converges 
to (Ti^/'c^ + {—iyRTL^i^)/\^ we get the conclusion of part (ii) for a subsequence first, then, 
by uniqueness of the hmit, for all L ^ oo. 

For part (iii), it suffices to show that there are no sequences (Ql, ip^) £ [uj—d^, u+dj^ x 
with ||'?/;^||j:^2 = 1 and ipL-^i^jiL), j = 0,1 in such that 

lim \\{HL-nL)ij''\\L^=0. (8.8) 

The spectral estimate: 

\\{Hl - nL)tp''\\L2 > dzst{nL,a{HL)\{no{L),ni{L)})\\tp'^\\L^ 

combined with ( I8.8P implies 

lim W^P^Wl^ = 0. (8.9) 

L— »oo 

In principle we can now employ the compactness argument in part (ii) to get 

lim W^j^Wh^ = (8.10) 

L— >oo 

which will contradict || ||/^2 = 1. More precisely, f l8.8p -( l8l9l) imply 

lim II (-A -uj-d, + Vl)iP^\\l2 = 

which, by repeating the argument after (18.11) with u replaced by + d*, gives 

lim W^'^ + TLi^M 

+ RTLip^+d^ 11^2 = 

L— »oo 

where ipuj+d* and ipuj+d* are eigenvectors of —A + vq corresponding to eigenvalue u; + rf^,. But 
the latter is not actually an eigenvalue, hence ipuj+d^ = and ipoj+d^, = 0. These show (18.101) 
and finishes the proof of part (iii). 

The proposition is now completely proven. 

Proposition 8.2 

fliooi + 2aioio - floooo < -7 < 0, and (8.11) 
" ^ as L t oo, (8.12) 



|aiooi + 2aioio — aooool^ 

These are now obvious from definition of aijki, the continuity of : x x t-^ 
and the convergence of ipj{L) to the translations and refiections of the single well 
eigenvector, see Proposition 18. II part (ii). 
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Figure 6: This figure shows the nonlocal analog of Figure [3l Panel (a) shows the dependence 
of the (absolute value of the) coefficients of the finite-dimensional approximation on the 
non- locality parameter e (e = denotes the "local" nonlinearity limit). The solid (blue) 
line denotes aoooo; the dashed (green) aim, the dash-dotted (red) aoiio; while the thick solid 
(black) one denotes aoioi- Panel (b) is analogous to panel (a) of Figure [3l but now shown 
for the non-local case, with the non-locality parameter e = 5. Finally, panel (c) shows 
the dependence of the critical point of the finite dimensional bifurcation (A'^,f2j,), on the 
non-locality parameter e. 
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